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Q\ . ABSTRACT 
Q\ ■ 

t— I ' We study the potential of the Pierre Auger detector for horizontal air showers 

(— i , initiated by ultra high energy neutrino. Assuming some simple trigger require- 

^ ' ments we obtain measurable event rates for neutrino fluxes from AGN, from 

, topological defects and from the interactions of cosmic rays with the microwave 

background. 

! 1. Introduction 

> ; 

It has been known for a long time that deeply penetrating high energy particles 
O ■ such as muons and neutrinos initiate horizontal air showers that can be detected at 

| : ground level B. Sinee the interaction length for muons and neutrinos in the atmo- 

sphere is larger than the whole atmospheric depth, they have roughly equal probability 
to interact at any point in the atmosphere. On the other hand, the rate of air showers 
due to the hadronic particles, that constitute the bulk of the cosmic rays, decreases 
very rapidly with zenith angle as the atmospheric depth rises from about 1000 g cm~ 2 
in the vertical direction to close to 36000 g crrT 2 horizontally. The electromagnetic 
component of air showers started by electrons, photons and hadrons gets absorbed 
well before reaching the Earth's surface and only the muon component of the shower 
survives for sufficiently large zenith angles. A detector that is able to identify the 
electromagnetic component of air showers is then capable of identifying horizontal 
showers induced by such penetrating particles. Such an array will mainly trigger on 
horizontal showers that initiate at the appropriate depth so that the shower is close 
to shower maximum when it reaches the array. 

The recent agreement between calculations of diffuse neutrino fluxes from Active 
Galactic Nuclei (AGN) has raised a lot of expectations for neutrino telescopes that 
are currently under development and construction. Although horizontal showers have 
ruled out an early prediction of neutrino fluxes from AGN these fluxes extend to 
the PeV region where the corresponding horizontal shower rate is very close to that 
expected from hard bremsstrahlung of the conventional atmospheric muon flux (pro- 
duced in it and k decays). Horizontal showers are currently being studied by several 
ground arrays because they should provide complementary information on prompt 
muon and neutrino production in the atmosphere which can be related to production 
of charm & and of cosmic ray composition around the knee it. It is accepted that the 



most appropriate technique for neutrino detection consists on detecting the Cerenkov 
light from muons or showers produced by the neutrino interactions in water or iceS*. 

The situation is however different for still higher energy neutrinos where the 
project to build two 3000 km 2 particle arrays one in each hemisphere (Pierre Auger 
Detector) may play an interesting role. The project is discussed in a separate article in 
these proceedings©. The reference design combines an array of particle detectors and 
an air fluorescence device similar to Fly's Eye to detect cosmic ray air showers of en- 
ergies above 10 19 eV. The proposed particle detectors are water Cerenkov tanks, very 
appropriate for detecting particles arriving horizontally. The detector will be most 
efficient for high zenith angle showers of energy above 10 19 eV when a large number 
of detectors of the array register significant signals. The electromagnetic component 
is separated from the muon component on the basis of the individual muon pulses 
that stand out of the average signal produced by the electromagnetic component of 
the shower. Neutrino predictions of such energies include those from interactions of 
the cosmic rays with the cosmic microwave background which have a solid foundation 
and would be of enormous value to establish the Greisen-Zatsepin-Kuz'min cutoff, as 
well as more speculative sources such as topological defects!! and primordial black 
holes!) . 

2. Large Showers from Neutrinos 

Neutrinos produce showers in most interactions with the atmosphere but the show- 
ers are of different nature depending on the process in consideration. In the interac- 
tions the target nucleons break up and the debris behaves as a group of hadrons that 
results in a shower similar to those induced by regular hadronic cosmic rays. Such 
shower is produced in both neutral and charged current interactions. If the neutrino 
is of electron flavor there will be an additional shower produced by the electron at 
the leptonic vertex of charged current interactions. This shower is of electromagnetic 
type, somewhat narrower and with a smaller muon content than the hadronic showers. 
The resulting shower for the charged current interactions of electron neutrinos is thus 
a superposition of two parallel showers one of hadronic type and the other electromag- 
netic. For high energy neutrino interactions the average fractional energy transfer to 
the nucleon in the lab frame (y) is < y >~ 0.2 so the electromagnetic shower carries 
on average 80% of the neutrino energy and is therefore most important. The interac- 
tions of the neutrinos with the electrons have in general much smaller cross sections 
and can be disregarded except for the resonant electron-antineutrino electron interac- 
tion which dominates just for neutrino energies around the resonant value of 6.4 PeV. 
In that case the character of the shower depends on the disintegration channel of the 
produced W boson in the s-channel. 

The showers regardless of their character can be detected by a particle array if 
they are initiated at an appropriate distance by a neutrino with sufficient energy. 



They will resemble ordinary air showers the main difference being that horizontal 
showers develop in a more uniform atmosphere. Electromagnetic showers should 
have lateral and depth distributions according to standard parametrizations when 
the corresponding lengths are measured in depth (g cm -2 ). Similarly we assume that 
hadronic showers will be not too different to ordinary cosmic ray showers. Large zenith 
angle showers of energies above the array threshold can be detected if a suitable trigger 
is selected, provided the plane of the array intersects the shower at a location where 
the number of particles is close to its maximum. This is a conservative statement 
since the particle detectors of the array are more closely distributed in the transverse 
plane to a near horizontal shower. The triggering to be developed will be however 
very different from the standard trigger for vertical showers, in particular the relative 
timing of the signals in adjacent detectors will reflect the shower propagation across 
the array. 

3. Event Rates 

For a neutrino flux d& u /dE u interacting through a process with differential cross 
section da/dy, where y is the fraction of the incident particle energy transferred to the 
target, the event rate for horizontal showers can be obtained by a simple convolution: 

$ sh [E sh > E t h] = N aPmr / dE sh / dy -r^{E v ) —{E v , y) A(\, S v ) (1) 

JE th Jo dE v dy 

where N a is Avogadro's number and p air is the air density. The energy integral 
corresponds to the shower energy E s h which is related to the primary neutrino energy 
E v in a different way depending on the interaction being considered. A is a geometric 
acceptance which contains the volume and solid angle integrals for different shower 
positions and orientations with respect to the array. 

3.1. Acceptance 

We define the effective acceptance A as the integral over volume and solid angle in 
dVl = d(pd(sin 6) where <fi is the azimuthal angle in the array plane. It depends on the 
energy transfer to the shower and on the type of shower produced in the interaction. 
We take the shower axis to go through the array and assume showers are large enough 
to trigger when they start at an adequate point. The effective area is then simply 
given by S sin 9 where S is the area covered by the array and 6 = 90° — 9 Z is the 
angle between the neutrino arrival direction and the array plane, (complementary to 
the zenith angle). To obtain the effective volume it must be multiplied by a "depth 
interval" . For a given neutrino direction and impact parameter the depth interval is 
basically the range of positions of the interaction point that will trigger the array. 
As a given shower is moved through all possible first interactions points it spans 
an infinite cylinder which we refer as a shower-tube. Such tube intersects the array 



plane in an ellipse with a major axis given by q = IrjsinQ where r is the radius 
of the shower-tube. We can take the projection of this length onto the shower axis 
as the depth interval, which is equivalent to demanding that the shower maximum 
intercepts the array. This is conservative because it ignores the shower length which 
increases the range of allowed positions for the first interaction point. Since the depth 
interval cannot exceed the length of the array, W, we take the minimum of q and the 
average length of the array W. 

We calculate A integrating this volume over the possible solid angle orientations 
of the shower, dQ = 2irdsin6, and restrict the integration to horizontal showers i.e. 
0° < 6 < e max ~ 20°. 

A = Sxen \ , [(sin£)sin^-,+ / \(smP) sinPW 
\ Jsin^o sin^ •>' 



S x 2n r(2sm8 max — sin8i) (2) 



For 6 < 9i = sin 1 (2r/W), the intersection of the tube reaches W, its maximum 
value. 



3.2. Sensitivity to Neutrino Fluxes 

It is now a matter of substituting reasonable values for the parameters to get 
an estimate for the acceptance. For S = 3000 km 2 the "diameter" of the array is 
approximately D = 65 km. For W we should take the average length across the array 
for all possible impact parameters, we obtain W 0.70D ~ 45 km. The acceptance 
scales with the tube radius which we take as r = 1.5 km, the separation between 
the individual detectors of the array. Ordinary cosmic ray showers are expected to 
give measurable signals in detectors that are this distance away from the shower axis. 
The showers we consider here have shower maximum intercepting the array plane 
so they should have similar particle densities. We can now obtain B\ = 0.07 rad 
and an acceptance of A = oof// /V which when multiplied by an air density 
p air ~ 1.1 10~ 3 g cm~ 3 gives 2 10 7 kT sr. Neutrino detectors in planning aim towards 
an active volume in the range of 1 km 3 0* . Their effective volume is enhanced because 
of the long range of the energetic muon produced, but for electron and tau neutrinos 
they have to collect the the Cerenkov light from the showers they produce. If their 
energy is well above the PeV region the Earth will be opaque to these neutrinos and 
the corresponding acceptance of a 1 km 3 detector for contained events is at most 
6 10 6 kT sr, illustrating how the Pierre Auger project may come into play. 

In order to obtain a rough estimate of the rate of horizontal showers above 10 19 eV 
produced by a given neutrino flux we can simply take the product of the neutrino flux 
above 10 19 eV, the total cross section a and the acceptance A. For charged current 
neutrino electron interactions all the neutrino energy is transferred to the shower. The 
cross section corresponding to charged current neutrino interactions at this energy is 



uncertain because of the unknown behavior of the structure functions at low x and 
high Q 2 which take part in the calculation. Extrapolations of the structure functions 
lead to cross sections in the a = 1.3 — 4 10~ 32 cm 2 ranged. If we take the extreme 
neutrino fluxes from topological defect models & divided by a factor of 2 to account 
for electron neutrinos, the integral neutrino flux ranges from $ = 10~ 16 [cm 2 s sr]^ 1 
for the model with p = 1.5 to $ = 4 10~ 14 [cm 2 s sr]~ l for the model with p = 0. 
We obtain 2 — 5 10~ 8 s _1 , in the range of one event per year for the lowest flux and 
about 400 times that for the highest. The result is extremely encouraging because the 
calculation is very conservative. There are several issues that will rise the event rate: 
we have ignored the neutral current interactions and muon neutrinos all together, 
the cross section and the acceptance integral should both increase with energy, large 
showers which are very horizontal may trigger the array even if their axis falls outside 
the array area and it is also possible that the particle arrays have a lower threshold 
for horizontal showers. 
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Fig. 1. Neutrino flux predictions in the EeV range. 

The neutrino flux predictions for topological defects have been normalized in a 
maximal way, assuming the observed highest energy cosmic ray spectrum is due to 
the topological defects themselves but it may be that such fluxes are close to ten 
orders of magnitude below @. It should be stressed that there are solid predictions 
for neutrinos produced in the cosmic ray interactions with the cosmic microwave 



background responsible for the GZK cutoff. In the range 10 19 — 10 20 eV they only 
differ from topological defects by less than one order of magnitude (see Fig. 1). It is 
well possible that such interesting events become accessible to the two detector arrays 
planned. 

3.3. Full Event Rate Calculation 

For more realistic calculations, which are in progress, triggering details become 
important. We estimate such effects demanding that a number of consecutive particle 
detectors in a row have an electron density above a fixed value. The integral for 
the effective acceptance can be calculated numerically using parametrizations of the 
lateral distribution functions for both electromagnetic and hadronic showers. The 
results obtained in both cases are quite similar. As the threshold electron density is 
decreased the detector increases its acceptance for showers of lower energy because 
these showers are small and have to be extremely well aligned with the detector 
rows in order to trigger. However for the larger energy showers the acceptance does 
not change much as the threshold is lowered. Preliminary results for three triggering 
requirements are illustrated in Fig. (2) reflecting the stability of the results for shower 
energies above 10 19 eV. 
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Fig. 2. Acceptance calculation for three trigger models 



To calculate event rates we use two sets of structure functions MRS(G) © and 
GRV0. For the first we extrapolate to low x beyond validity of the parametrization 
using the slope of xq(x), where q(x) is the standard parton distribution. The second 
set, GRV, can be cautiously used on its own for low x. Fig. (3) shows both predictions 
in comparison with a data point obtained from HI collaboration in HERA El. We 
take three neutrino fluxes for reference calculation. We use the lowest prediction of 
ref. 0!) for neutrinos produced in the decay of topological defects with p = 1.5. This 
flux is very flat and would dominate the neutrino sky for energies above E v = 10 17 eV. 
We take the upper limit of the band calculated in ref. and the prediction of the 
neutrino fluxes from cosmic ray interactions with the cosmic microwave background 
calculated in ref. integrated up to redshift z = 2. Fig. (1) illustrates these fluxes 
compared to other predictions setting the scale of the sensitivity of the Pierre Auger 
project to high energy neutrino fluxes. We approximate the electron neutrino flux 
to be a factor of two below the muon neutrino for all three cases. This ratio can be 
naively expected from the number of channels in the decays of pions. The results are 
shown in table 1. 




Fig. 3. Comparison of two neutrino cross section predictions in the EeV range 



The acceptance curve shown in Fig. (2) is a continuous function of shower energy 
and when it is combined with the AGN flux prediction it can give measurable rates. 
Because these fluxes are typically of PeV energies they will produce small showers 
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0.9 




TD p = 


1.5 


1 


26 


9 


0.1 


51 


17 



Table 1: Yearly neutrino event rates for diffuse fluxes from AGN, for neutrinos from 
the interactions with the cosmic microwave background (CMB) and for topological 
defects in the model described in the text (TD). 

compared to the typical showers detected in the Auger detector and the showers will 
have to run well aligned with a row of particle detectors to trigger. These showers 
will undoubtedly produce signals that are very different from those of typical showers. 
The low energy part of these curves is very sensitive to the triggering conditions and 
there are large differences between event rates for different trigger models. This is 
not the case for the topological defect fluxes and for the flux from interactions of 
cosmic rays with the cosmic microwave background. These fluxes are much flatter 
and hence the horizontal shower rate peaks for neutrino energies in the region where 
the acceptance integral is fairly stable strengthening the results obtained. 

4. Conclusions 

The Pierre Auger project can be made sensitive to ultra high energy neutrino 
fluxes through horizontal showers if an appropriate trigger is implemented. Its accep- 
tance for detecting contained neutrinos events of energy above E u ~ 10 19 eV will be 
of the order of other neutrino telescopes in planning. The peak of horizontal shower 
acceptance for the Pierre Auger Project is at energies about 10 19 eV, a lot higher 
than the optimal region for AGN neutrino detection, for which the conventional ap- 
proach to detect neutrinos is best suited. The Pierre Auger Project is best suited 
for detection of neutrinos from interactions of cosmic rays with the cosmic microwave 
background and from the decay of topological defects. The event rates expected under 
some simplifying models for the trigger are high enough to be observed. 
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